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Abstract 

A Mathcmatica package for color summed calculations in QCD (SU(iV c )) is presented. Color contractions 
of any color amplitude appearing in QCD may be performed, and the package uses a syntax which is very 
similar to how color structure is written on paper. It also supports the definition of color vectors and bases, 
and special functions such as scalar products are defined for such color tensors. 

1. Introduction 

; 

With the LHC follows an increased demand of exact calculations in QCD involving many color charged 
Q_[ partons. Due to the non-abelian nature of QCD this poses a nontrivial computational problem. 
Q ■ This is the issue the Mathematica® package ColorMath is set out to tackle. The main feature of 

ColorMath is thus the ability to automatically perform color summed calculations starting from a QCD 
color structure which is expressed using a syntax very similar to how the color structure would have been 
written on paper. 

ColorMath allows both partial and full dummy index contractions, and thus deals with color structures 
having an arbitrary set of free indices. In this sense it is thus more general than on the other hand it 
is a pure SU(7V C ) tool, rather than a high energy physics general purpose tool, such as @, ColorMath 
works for an arbitrary N c > 2, and for arbitrary trace convention of the squared SU(7V C ) generators, Tr. It 
may be used by just giving the color structure in the appropriate form and afterwards running CSimplify 
to contract all repeated indices. For example, consider q\ q 2 ~^ 93 94 via gluon exchange in the t- and u- 
channels. Writing down the amplitude as 

Amp = St^\ 2 t^\ 3 +Tt^ 1 q3 t^ i q2 (1) 

where S and T represent some s- and t-channel kinematics, and t^ ql q 2 the SU(iV c ) generator in the 
fundamental representation, otherwise typically denoted (t 9 ) ql q 2, we may calculate the squared amplitude 
(defined as in the scalar product eq. (O) using 

CSimplify[Conjugate[Amp] ReplaceDummyIndices[Amp]] , (2) 

immediately giving the answer 

HV 2 _ 1) 

T R]\T — ^ NcS - T ) S * + ( N c T - S)T*}. (3) 

Each squared amplitude can in principle can be calculated like this. However, ColorMath also facilitates 
the usage of color vectors and bases and has special functions for the calculations of scalar products and 
gluon exchanges. 

This paper, which is intended to be the main reference, is organized as follows: First a general intro- 
duction to the basic color building blocks is given in section [2] In section [3] the computational strategy is 



presented, whereas basic examples are given in section 01 and the usage of vectors and their functions are 
presented in section [3] In section [5] some remarks concerning validation and scalability are made and in 
section [7] conclusions are drawn. 



2. QCD building blocks 

From a color perspective, the QCD Lagrangian is built out of 



jo a 
i K 3 



quark-gluon vertices, - — * — = (t a ) l j, (4) 

and 



a, . 

ibc 



triple-gluon vertices, ^^TTo^ c = C ' (5) 

where we follow the convention of reading of the fully antisymmetric structure constant indices in counter- 
clockwise order. The four gluon vertices can be rewritten in terms of (one gluon) contracted triple-gluon 
vertices, and thus need no special treatment. 

Due to confinement, we never observe individual colors and it therefore suffices to calculate color 
summed/averaged quantities for making predictions in quantum chromodynamics. We may thus constrain 
ourselves to treat QCD amplitudes carrying a set of external indices with values that need never be specified, 
as they are always summed over in the end. 

In principle, the SU(iV c ) generators, along with a delta function for indicating a quark and an anti- 
quark color singlet, S ql q 2, constitute a minimal set of objects for treating the color structure in QCE0. 
For convenience, and for performance reasons, it is, however, useful to define a larger set of objects. The 
complete set of basic building blocks for carrying color structure used by the ColorMath package is given in 
tabled] 

Apart from a delta function in quarks indices, S ql q 2, a delta function in gluon indices, denoted by A^ sl ' fl2 ^, 
is also defined. Note the Mathematica List brackets in {gl,g2}. The gluon delta function can alternatively 
be expressed using 

5 W0«/ = Tr * tr ^ 92 ] = T « 59192 * ° {9h92i = TRA{9h92} ® 

where Tr (TR) comes from the normalization of the SU(iV c ) generators, and is typically taken to be 1/2 
(the Gell-Mann normalization) or 1. A rescaling of the normalization of the SU(iV c ) generators can always 
be absorbed into a normalization of the strong coupling constant. To allow for arbitrary normalization Tr 
is kept as a free parameter, denoted by TR, and may be defined by the user. Also the number of colors, 
denoted Nc in ColorMath, may be set by the user; by default both Nc and TR are kept as free parameters. 
Note the ColorMath notation for a trace over two gluons o^ 9l,g2 ^ in eq. ([5]). Similarly, a general trace over 
k gluons gl,g2,..gk is denoted o{ gl > g2 >-> gk } : and may be thought of as a closed quark-line with k gluons 
attached. 

The totally antisymmetric structure constants, which - along with and extra i — defines the triple 
gluon vertices, are denoted ffs 1 ^ 2 ^ 3 ). Similarly, the totally symmetric "structure constants" are defined as 
di gl ' g2,g3 ^ . Recall that, starting from the commutation (anticommutation) relations 



[t gl ,t g2 ] = if glg2g3 t g3 , {t 9 \t 92 } = d 9l92g3 t 93 , (7) 



1 To enhance the similarity with usage inside Mathematica, we here use the somewhat unorthodox notation ql etc. for single 
quark and gluon indices. 
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Table 1: Along with the number of colors, Nc, and the trace of an SU(iV c ) generator squared, TR, the basic building blocks arc 
as below. Note that o^ 1 ^ 2 - ••■•9 fe } represents a trace over gluons with indices gl...gk, tr[t 9l t 92 ...t 9fe ], and that t* 91 ' 92 '— ' gk ^ 1 q 2 
represent the ql, q2-component in a trace over generators that has been cut open, (t gl t 32 ...t 9k ) ql q 2- For convenience, also the 
totally symmetric "structure constants" (j{ 9l > 92 > 93 } are defined. Note the Mathematica FullForm. ColorMath is built on 
pattern matching, and it is therefore essential that the expressions have the correct FullForm. In particular, Power may 
not be used instead of Superscript. To get the right FullForm it is recommended to use the function form, which is just 
a function returning the corresponding ColorMath object. 



Pictorial representa- 
tion 



ColorMath 



Function form 



Mathematica FullForm 



ql q2 
ffl 52 

gl gl gk 
?1 11' I ?2 



gl g2 



gk 



c 



s ql q2 

&{ghg2} 

f{gl,g2,g3} 
d {gl,g2,g3} 

t{gi,— ,gk}qi 
{gi,-,gk] 



q2 



<5[ql,q2] 
A[gl,g2] 
f[gl,g2,g3] 

d[gl,g2,g3] 

t[{gl,...,gk},ql,q2] 
o[{gl,. . . ,gk}] 



Subscript [Superscript [\ [Delta], ql], q2] 
Superscript [\[CapitalDelta] , List[gl, g2]] 
Superscript [f, List[gl, g2, g3]] 

Superscript [d ,List[gl, g2, g3]] 



Subscript [Superscript [Superscript [t , 

List[gl,...,gk]],ql],q2] 



Superscript [o, List[gl, . . . , gk]] 



the structure constants can be rewritten in terms of traces over SU(iV c ) generators, 

if 91 92 93 /d gl 92 93 = —\tr{t 9l t 92 t 93 )Ttr{t 92 t 9l t 93 )] 
Tr 

= [(t 9l ) ql q 2(t 92 y 2 q3 (t 93 ) 93 ql t (t s2 ) ql q 2(t sl ) q2 q3 (t 93 ) 93 ql ] 



(8) 



In ColorMath notation t is used to denote an open quark-line, and the above expression is written similarly, 



Pictorially this represents 



53} _ {gl,g2,g3} {g2,gl ig 3}\ 

TR 



(9) 



J-Ct^ 1 ^ 1 n t^ 2 !"? 2 q t {93}g3 ,{g2} q l ,{ 5 l}q2 .{ 9 3} q 3 ) 

TR^ ?2 l 9 3 1 ql) 




1 






1 

Tr 




(10) 



where i is included in the triple gluon vertex. The rationale for putting the gluon index in the SU(iV c ) 
generators t^ g2 ^ ql q 2 inside a Mathematica List in eq. ^ is to allow for the natural extension of having 
many gluons attached to an open quark-line, thus 
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The left hand side has several advantages compared to the right hand side. An open quark-line with an 
arbitrary number of gluon indices can be written in a compact form with no dummy indices. This is not only 
more human readable, but also superior from a computational point of view, as it avoids the contraction 
of unnecessary dummy indices. A third advantage with the above notation is its direct correspondence to 
the trace type bases [4-ll2j. A basis (or spanning set) for the color space for a fixed set of external quarks, 
anti-quarks and gluons can always be taken to be a sum of products of open and closed quark-lines. 

In this context we also remark that the color tensors defined in table [T] arc color scalars, i.e., they are 
invariant under S\J(N C ) transformations. This imposes no restriction for our purposes as, for any QCD 
amplitude, the overall color structure, including both incoming and outgoing particles, always is a color 
singlet. As the basic building blocks are invariant, each tensor built out of these objects, i.e., each tensor 
needed for color summed calculations in perturbative QCD is a actually a color scalar. 

The scalar product on this vector space is given by summing over all external color indices, i.e. 

( Cl |c 2 )= cp a2 -c? a2 - (12) 

with a,; = 1, N c if parton i is a quark or anti-quark and ai = 1, N£ — 1 if parton i is a gluon. As long 
as the color structures in table [1] are multiplied by real coefficients the scalar product is actually real, which 
is easy to prove using the computational rules in the next section. 



3. Basic computational strategy 



Having defined all the color carrying objects, we turn to describing the basic strategy for carrying out 
calculations. 

Again we note that we need not treat the four gluon vertex as this can be rewritten in terms of three 
gluon vertices. To treat an arbitrary color structure in QCD we may always compute any squared amplitude 
in the following way: 

(i) Rewrite the triple gluon vertices using eq. ([5]). This results in a color structure which is a sum of 
products of open and closed quark-lines, connected to each other via repeated indices. 

(ii) Contract all internal gluon indices using the Fierz or completeness relation 
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(13) 



(iii) If present, remove Kronecker deltas (for quarks and gluons), and internal quark indices using 

gl grn gz gl g 2 gk gl gz 



ql 



9\ <?2 


gk 


q o 
q o 
q o 
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?2 
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(14) 



(iv) Use the tracelessness of the SU(iV c ) generators, and contract quark and gluon delta functions with 
repeated indices 




= Nl-l. 



(15) 



After this, all dummy indices have been contracted. A squared amplitude is thus represented by a rational 
function in Tr and N c , whereas a general amplitude is represented as a sum of products of closed and open 



{}?i 



<l2 



5 ql q2 , and {Bi,92,~,gk} ^ which 



quark-lines, i.e., as sums of products of t^ 9l ' 92, - ' gk ^ ql q 2, including t 
for two gluons may be rewritten as o^ 91 ' 92 } = TR A^ 91 ' 92 }. 

By applying these rules, any amplitude square and any interference term appearing in QCD can be 
calculated [a. Ill|. To successfully square an arbitrary QCD amplitude using the above set of rules we see 
that steps (i-ii) have to be performed while keeping the relative order. We also note that these rules increase 
the number of terms, whereas the rules (iii-iv) decrease the number of terms or keep it fixed. In order not 
to unnecessarily inflate an expression it may therefore be useful to apply the latter rules at any time during 
the computation. ColorMath utilizes this and tries to contract indices using (iii-iv) at any time, while rules 
(i-ii) are used only when needed, i.e., when the non-expanding rules fail. 

On top of the rules (iii-iv) there are other cases in which the contraction of a gluon results in at most 
one term. This will happen if 

(v) Two neighboring gluons attached to the same (closed or open) quark-line are contracted. In this case 
the result is simply cp — Tr(N 2 — \)/N c times the color structure where the two involved gluons have 
been contracted, 



gl gk gl 




gm 



N 2 



1 



gk gm gz 



c5 o 


o 


o 


Si q 


o 


o 


a - o 


o 


- o 


a o 


o 


o 


q o 


o 


o 



(16) 



(vi) Similarly, it is easy to show, using the Fierz identity, eq. (|13[) . that the contraction of two next to 
neighboring gluons results in only one term 




gm gn gz 
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(17) 



While the above rules are sufficient for squaring color amplitudes, it is sometimes advantageous to con- 
tract indices directly using the structure constants f gl g2 93 and d glg2g3 . Therefore rule (i) is not always 
automatically used for simplifying expressions. Instead, in addition to the above set of rules, a limited set of 
rules for contraction of repeated gluon indices occurring in the standard and the symmetric structure con- 
stants are implemented. More specifically, all rules for contracting two gluons in products of two symmetric 
or antisymmetric structure constants, for example, 

jgl il i2p2 il 12 = Tr ggl g2 ^ 

and all rules for contracting three gluons in products of three symmetric or antisymmetric structure con- 
stants, such as 



N 2 — 4 

^gl il i2^g2 il i3 jg3 12 13 _ c jglg2g3 

are implemented. For more intricate gluon contractions involving structure constants ColorMath can con- 
tract indices by first applying rule (i). 



4. Basic calculations 



In this section, basic functionality of ColorMath is explored. As always with Mathematica, the package 
has to be loaded before it can be used. With 

Get["/full/path/to/ColorMathl.0.m"] (20) 

where the version number is adjusted to the version in question, this can be done. To start using the package, 
it is, however, recommended to start from (a version of) the tutorial "ColorMathTutorial.nb" and modify 
its content according to the desired usage. 

The color contractions corresponding to the basic manipulations from section [3] are carried out using 
Mathematica Rules, i.e., a set of replacement rules based on pattern matching. As always, the rules may 
be applied using "Expr/.TheRules" , and may be applied repeatedly using "Expr//.TheRules" . 

The rules described in (iii-iv) are, along with eq. contained in a set of rules called SimpleRules, as 
they keep the expression at least as simple (in a term counting sense) as it initially was. Applying these to 
an expression thus tend to simplify it, for example 

t {al}ql q2 t {92}9 V A* 51 ' 52 * /.SimpleRules (21) 

results in TR(A< 9l ' ff2 >) 2 . 

To fully utilize the rules we (may) need to apply them repeatedly, 

t {9i}<H q2 t {92}q2 qi A {fli, ff 2} //.SimpleRules (22) 

giving the fully contracted expression (— 1 + Nc 2 )TR. 

The rules defined in (v-vi), acting on t and o, are similarly contained in OTSimpleRules, and the union 
of SimpleRules and OTSimpleRules, and with a few rules for rewriting closed quark-lines with zero to two 
gluons, are contained in AllSimpleRules. 

The special rules for gluon index contraction in structure constants, exemplified in equation (|18M19[) arc 
defined in FDRules. The complete set of rules are stated in table IA.2I 

Rather than thinking about how individual rules have to be applied, it is convenient to have a standard 
procedure for contracting color indices. This is embodied in the function CSimplify, which does what its 
name suggests; simplifies the color structure as far as possible. For color structures which does not contain 
structure constants this always implies contracting all repeated indices. 

For color structure involving the structure constants CSimplify first attempts simplification using the 
FDRules. If, after this, the expression still contains structure constants, the structure constants are by 
default rewritten in terms of traces using eq. ©, and the indices are fully contracted, resulting in a sum 
of products of open and closed quark-lines. Sometimes it may, however, be desirable not to rewrite the 
structure constants, as expressions may be more compact if they are kept. This can be achieved by using 
the option RemoveFD— > False. 

The most useful set of functions are given in table IA.3I Apart from CSimplify we especially note the 
function ReplaceDummylndices for replacing all repeated indices in a color structure with a new unique 
set of color indices. 

For the purpose of calculating amplitudes square, we need, as in eq. (fT2")> . the complex conjugated version 
of the color structures. We note that 

Conjugate[5 9l 92 ] = 5 q2 qX 
Conjugate^ 1 - -^ 1 q2 ] = t ^ -^ 2 gl 

Conjugate^ 1 '-' 9 **] = o^ fc '-' ffl > (23) 

whereas A< sl ' 92 > , flsLs 2 ^ 3 ! and d {gl > 92 ' 93} are real. This is implemented in ColorMath via redefinition of 
the Mathematica's built in function Conjugate. With this in mind, we are ready to perform calculations of 
the type in eq. J2]). 
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Finally we remark that to Mathematica A^ 1 '«' 2 > ^ A^ 2 ^ 1 }, and similarly, o^ 2 ^ 3 -^ 1 } ^ {<?i,s2, 9 3, 9 4} 
etc. For this reason a function Sortlndices, which writes indices in Mathematica default order is defined, 
s.t. for example 



SortIndices[o^ 2 ' ff3 ' 94 ' 9l >] = o^ 1 ^ 2 ^ 3 '^ . (24) 

This function is used by CSimplify, and by some of the rules in table [A~2"l 

Sometimes more detailed control over the calculation may be desired. For this purpose, and for internal 
usage, functions manipulating indices and probing the color structure are given in table [AT4l in | Appendix A| 

While each squared amplitude in principle can be calculated like this, ColorMath also offers more efficient 
tools for dealing with vectors in color space. 



5. Defining and using vectors 

For the purpose of studying color space it is often convenient to define a basis for the color space, and 
sometimes also projection operators. Both of these are examples of color (singlet) tensors, and ColorMath 
has a set of tools for working directly with such tensors. 

As an example, let us consider the color structure for q{q 2 — > <7394.<?e0- A basis for the color space may 
be written as [13[ 

Vectorl81 {gl _ ig2 _ >g3 _ )54 . lS 5_} := S ql q2 t {95}q4 q3 
Vector818 {9l _ )92 _ ig3 _ ;94 _ jff5 _ } := t^i\ 2 5"\ z 

Vector888s { g 1 _,g 2 _ i<?3 _ 1 g4_ )S 5_} := Module[{il, i2}, t^ ql q2 t^ q4 q3 d* 11 ' 12 '^] 
Vector888a {gl _, g2 _ i(? 3_ ig4 _ )S5 _ } := Module[{il, i2}, t* 11 ^ t { i2 >« 4 g3 If* 1 ' 12 '^] (25) 

where the basis vectors are labeled using first the overall multiplct of qiq 2 , then the overall multiplct of q^q^ 
and finally the overall multiplct of q^q^gz , which - due to color conservation - must equal the multiplet of 
qiq 2 , (implying that the notation is somewhat redundant). 

From a Mathematica perspective we note a few things. First, on the left hand side, we see that the 
indices inside the List in the Subscript are followed by underscore to indicate pattern matching. This is 
standard in Mathematica and makes it possible to use any symbol to denote the indices in later calculations. 
Then we note that the last two tensors in eq. ((25]) are defined using Module. This is to ensure that each 
time the tensor is used, it comes with a fresh set of dummy indices. This is also the reason why set delayed 
":=" is used. 

This basis is orthogonal since at least one set of partons transform under different representations in the 
various tensors. It is, however, not normalized. Finding the normalization is easy using ColorMath. To 
calculate for example Vectorl81 squared we could enter 

CSimplify[Conjugate[Vectorl81 {ql q2 ,q3,q4,g5}] Vectorl81 {ql,q2,q3,q4,g5}]> ( 26 ) 

but it is yet much easier to use the tensor functions for calculating scalar products. Instead we could simply 
write 



CDot [Vector 181, Vector 181] (27) 

resulting in Nc (—1 + Nc 2 ) TR. Having a basis we naturally want to calculate all norms (square). ColorMath 
has special functions for this as well. If we define 



2 In this user guide indices representing incoming quarks and outgoing anti-quarks are placed upstairs, whereas outgoing 
quarks and incoming anti-quarks are placed downstairs. Note, however, that we could as well have used the opposite convention. 



OurBasis = { Vector 181, Vector818, Vector888s, Vector888a} 

we may calculate the squares of the basis vectors using 



(28) 



CDot [OurBasis]. (29) 

This results in a List containing the scalar products between each basis vector and itself 

, 2 . , 2 , 2 (4 - 5Nc 2 + Nc 4 ) TR 3 , „. „1 

Nc (-1 + Nc 2 ) TR, Nc (-1 + Nc 2 ) TR, — i — '- , 2Nc (-1 + Nc 2 ) TR 3 \ . (30) 

Often non-orthogonal bases are used, in particular, this is the case for the trace type bases. In this 
case, to square a general amplitude, all the scalar products between all the basis vectors are needed. For 
calculating the scalar product matrix the command CDotMatrix may be used. This returns a List of Lists, 
i.e. a matrix, where the ij-th element is the scalar product between the i-th and j-th vector in the list 
of (basis) vectors. For larger vector spaces with complicated scalar products, it may be desirable to get 
progress information on the calculations. This can be obtained by setting the option Verbose — > True for 
CDotMatrix, 

CDotMatrix[OurBasis, Verbose -)• True] . (31) 

Similarly, it may be useful to be able so simplify potential (normalization) roots assuming a large Nc. This 
can be done by using the option NcMin. The scalar product functions, along with their options, are listed 
in table IA.5I 

Often it is also of interest to investigate the effect of gluon exchange on the color structure expressed in 
a basis, i.e. starting in a basis vector j, what is the effect on the basis vector of exchanging a gluon between 
partonl and parton2. This is useful both for soft gluon resummation and for one- loop corrections via gluon 
exchange. The result can be expressed in terms of a matrix whose clement ij is the i-th component resulting 
after such an exchange in the initial vector j. This is calculated by the function 

CGamma[Basis, partonl, parton2] (32) 

where Basis is a List of basis vectors, defined using the syntax in eq. (|25p and partonl and parton2 are the 
numbers of the partons in the basis vectors, i.e., in this case numbers between 1 and 5. The sign conventions 
are such that quark-gluon vertex always comes without additional signs, and the triple-gluon vertex have the 
indices appearing in the order: external index, internal dummy index and index of the exchanged gluorjf]. 
For example, we may calculate the effect of gluon exchange between parton 1 and 3 in the above basis using 

CGamma[OurBasis, 1, 3, Verbose — > False, BasisType — > OrthogonalBasis] . (33) 

Here we have supplied optional information about the basis type, that the basis is orthogonal, to speed up 
the calculations. By default CGamma docs a few consistency checks. It checks that the vector resulting 
after gluon exchange, when squared, has the same value as the basis decomposed version. For orthonormal 
bases (BasisType — > OrthonormalBasis) it is also checked that the resulting matrix is symmetric 
These checks may, however be turned off (MakeChecks — > False). The set of options, with default values 
are listed in table IA.5I 

6. Validation and scalability 

The computational rules and functions in this package have been used for calculating the three gluon 
projection operators presented in [TEj] . This imposes highly nontrivial consistency checks, as it is verified 

3 The sign conventions thus differ from the typical eikonal choice. 
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that every projector square equals itself, which at intermediate steps often involve several ten thousand 
terms. Additional consistency checks on the color contraction rules have been made by using the CGamma 
function which checks that vectors square and basis decomposed vectors square agree, and by changing the 

order in which the color structure contraction rules are applied. The scalar product matrices have also been 

i — i 

compared to the ColorFull code 16| for tree level trace bases with up to six partons out of which one parton 
is a quark and one an anti-quark. Selected results have been compared against 13J, and the package has 
been tested both in Mathematica 7 and 8. 

The computational effort needed for exact treatment of the color space grows very quickly with the 
number of partons. The dimension of the vector space grows roughly as a factorial in the number of gluons 
plus qg-pairs [TBI (strictly speaking an exponential for finite N c in a multiplet basis). The computational 
effort for ColorMath, or any program operating by direct manipulation of quark-lines, tend to grow roughly 
as the square of this, as the bases are non-orthogonal. ColorMath (in its current form) is thus rather intended 
to be an easy to use package for calculations of low and intermediate complexity than a competitive tool for 
processes with very many colored partons. 



7. Conclusion 

In this paper a Mathematica package ColorMath for performing color summed calculations in QCD is 
presented. This package allows for simple evaluation of QCD color amplitudes which are expressed in a 
format which very much resembles how the color structure would have been written on paper, see table [TJ 
The idea is that the user - for simple cases - just should give the expression, and then run CSimplify[Expr] 
rather than Simplify [Expr]. The package is based on advanced pattern matching rules, and a list of rules 
is given in table IA.21 whereas functions acting on color structures are given in table IA.3I 

For calculations of intermediate or high complexity it is often beneficial to use a basis for performing 
color space calculations. ColorMath allows for definition of color tensors of form Cl{ il _ i2 _ ik_\ '■= ■ ■ ■, 
carrying an arbitrary set of quark, anti-quark and gluon indices. Special functions for calculating scalar 
products, and investigating the effect of gluon exchange, are given in table IA.5I ColorMath is, however, 
not intended for hi gh speed calculations involving many colored partons. For this purpose a separate CH — h 
package is written [la ]. 
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Appendix A. Rules and functions 



In this appendix rules and functions for index and color structure probing are stated. First the basic 
rules for manipulations are stated in table IA.21 then the most important functions, operating directly on 
the basic color objects from table Q] are listed in table IA.31 whereas special functions for indices and color 
structure probing are stated in table IA.41 and the functions dealing with color vector objects are stated in 
table IA.5I 



Table A. 2: The complete set of rules used by ColorMath to contract indices. These rules are applied as always in Mathematica, 
by using "Expr/.TheRule" . To ensure that all indices which can be contracted by the rules actually are contracted, apply 
the rules repeatedly, "Expr//.TheRule" . Note that all rules have names ending with Rules. 



Rule 


Effect 


AllSimpleRules 


AllSimpleRules is the set of all rules involving A^ 9l < 92 >, 8 ql q2 , o^ 1 '-' afc > and 
f.{gi,...,gk}qi^ 2 w hich do not increase the number of terms, i.e., the union of 
SimpleRules and OTSimpleRules. 


ExpandThenRules 


The rules ExpandThenRules first uses Expand, then applies OTThenAllSim- 
pleRules and finally restores default index order with the Sortlndices function. 


FDRules 


Rules for gluon contraction for terms involving up to three structure constants, 
f{gi,g2,g3} Qr d { fl i )fl 2, fl 3} ) ; e ^ expressions of form eq. (JTSJl and eq. (US]). The 
terms resulting after contraction are expressed in terms of structure constants 
and gluon deltas, A* 91 ' 92 *. 


FDToORules 


Rules for replacing structure constants, fi 9l >f 2 > 93 l and d^ 91,92 ' 93 ^, with sums of 
closed quark-lines o'f 91 ' 92 ' 93 '' using the intermediate expression in eq. 


FDToTRules 


Replaces structure constants, fl 9l > 92 ' 93 J and d^ 9l,92:93 ^, with a sum of products 
of SU(7V C ) generators t^ gl ^ ql q 2 using the last form in eq. (j9]). 


OTGluonRules 


Rules for contracting repeated gluon indices in ois 1 '—^} and t^ 9l, '"' Bk ^ ql q 2- 
The Fierz identity, eq. (1131). is included in these rules. 


OTSimpleRules 


Rules for contracting neighboring and next to neighboring gluons in closed and 
open quark-lines, {9h-,gk} an( j t {gi,...,gk} q i eq ^ and eq jy^ 


OTThenAllSimpleRules 


The rules OTThenAllSimpleRules first applies OTGluonRules and then re- 
peatedly AllSimpleRules. 


OTToTRules 


Rules for replacing open and closed quark-lines, o^ sl '""' sfe ^ and t^ 9l, "' gk ^ ql q 2, 
with products of SU(iV c ) generators t^ 9l ^ ql q 2- 


RemoveOORules 


Rule for simplifying closed quark-lines with gluons, = Nc. 


RemoveOTolORules 


Rules for simplifying closed quark-lines with or 1 gluons, = Nc, o^ 9 ^ = 0. 


Remove0To2ORules 


Rules for simplifying closed quark-lines with to 2 gluons, = Nc, o^ 91 ^ = 
and o* 91 ' 92 } = TRA* 91 ' 92 *. 


SimpleRules 


Basic rules for quark and gluon contraction. These rules involve A^ 9l ' ff2 ^, 8 ql q 2 
or quark contraction, and never increase the number of terms. These rules thus 
contain the rules in (iii-iv) as well as in eq. EJ). 
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Table A. 3: The most useful functions, to be used on any expression carrying color structure of the form given in table ^ 
Recall to relabel dummy indices, manually, or by using ReplaceDummylndices, and to place quark indices that are to be 
contracted such that one index sits upstairs and one downstairs. Furthermore, check that the correct FullForm is used, cf. 
tabled The function ^Vhatls Wrong may be useful for identifying common input mistakes in the basic color objects. 



Function with usage 


Effect 


CSimplify[Expr] 

Options with default value: 

RemoveFD — > False 


The most general function for simplifying color structure. If Expr 
contains structure constants, ffs 1 ^ 2 ^ 3 } or d {gl ' g2 > g3} , FDRules 
are hrst applied. If, after this, the expression still contains struc- 
ture constants they are - be default - removed using FDToORules, 
and all repeated indices are subsequently contracted by repeat- 
edly using AllSimpleRules, then OTThenAllSimpleRules, and 
finally ExpandThenRules. For color structure containing struc- 
ture constants, it could happen that it is desired not to replace the 
structure constants. This can be achieved by setting the option 
RemoveFD — ^ False. 


GluonContract[Expr, Gs] 


Contracts a set (List) of gluons Gs = {gl,...,gn} or a single 
gluon Gs = gl, in the expression Expr, while leaving other indices 
uncontracted. This function is intended for quark-lines and will 
replace structure constants with quark-lines. 


ReplaceDummylndices [Expr] 


Replaces the dummy indices in Expr with a new set of unique 
dummy indices. 


Sortlndices [Expr] 


Sorts the gluon indices appearing in A^ 1 ' 9 ^, fis 1 '^ 2 'f 3 I or 
d { ff i,s2, S 3} and Q { g i,...,gk} such that they stand m Mathematica 

default order. This is needed to ensure that one color structure 
only is represented in one form. 


SplitConstAndColor[Expr] 


Splits an expression Expr into a List of Lists of color struc- 
tures and corresponding multiplicative factors {{constants 1, color 
structure 1}, {constants 2, color structure 2},...}. This is done by 
first expanding the expression and then splitting the terms sepa- 
rately. 


WhatIsWrong[Expr] 


Checks if anything is obviously wrong with an expression, for ex- 
ample if Power is used instead of Superscript or if gluon indices 
are placed downstairs. The check is performed by first expanding 
the expression, and then checking each term. Read error messages 
from above. 
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Table A. 4: The below set of functions are used for probing the color structure, and identifying indices. Concerning the 
convention of what (lower or upper) quark-indices are referred to as incoming quarks and outgoing anti-quarks, as opposed 
to incoming anti-quarks and outgoing quarks, ColorMath does not specify any preference, but simply refers to the indices as 
upper quarks and lower quarks when needed. 



Function 


Effect 


Alllndices [Expr] 


Returns a List of all (external and dummy) indices in Expr. 


ContainsColor [Expr] 


Returns True if the expression Expr contains any color structure, i.e., 
any of the terms in table [TJ and False otherwise. 


ContainsFD[Expr] 


Returns True if the expression Expr contains structure constants and 
False otherwise. 


ContainsGluonDelta[Expr] 


Returns True if the expression Expr contains a gluon delta-function, 
A^ g1 ' 92 ^ , and False otherwise. 


ContainsO [Expr] 


Returns True if the expression Expr contains closed quark-lines, 
{gi,—,gk} ^ an( j p a j se otherwise. 


ContainsQuarkDelta [Expr] 


Returns True if the expression Expr contains a quark delta-function, 
S ql q 2, and False otherwise. 


ContainsT [Expr] 


Returns True if the expression Expr contains open quark-lines, 
t { g i,...,gfc} g i^ and Falge otherwise. 


Dummylndices [Expr] 


Finds the dummy indices in an expression Expr, by first expanding it 
and then finding all dummy indices in all terms. 


Gluonlndices [Expr] 


Returns a List of all gluon indices (external and dummy) in the expres- 
sion Expr. 


Lower Quarklndices [Expr] 


Returns a List of all (external and dummy) quark-type indices placed 
downstairs, i.e., q2 in t {gl '-' gk ^ <l1 q2 and 5 ql q2 . 


UpperQuarklndices [Expr] 


Returns a List of all (external and dummy) quark-type indices placed 
upstairs, i.e., ql in t^- 9l '-' 9k ^ ql q2 and S ql q2 . 
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Table A. 5: Special functions for color structures expressed in the form of eq. d 2 5 D . i.e., the corresponding tensors need to 
be defined using the pattern matching underscores, potential dummy indices should be hidden inside Modules, and the 
evaluation should be delayed for expressions containing dummy indices, i.e. ":=" should be used. The vector indices should sit 
inside a List in the Subscript. Additional options may be supplied to these functions to control for example the verbosity. 
Thus, CNorm may for example be called using CDot[{Cl,. . . , Cn}, Verbose — > True], to get progress information or 
simply as CDot[{Cl,. . . ,Cn}] to use the default options. 



Vector function with usage 


Effect 


CDot[Cl, C2, Options] 
Options with default Vctluci 
NcMin -> 3 


Calculates the scalar product between two color ten- 
sors CI and C2. Before using CDot the color tensors 
C1.,i_.,j_ „,_; and C2 {ll _ ii2 _,...,ifc_} should have been de- 
fined. For simplifications (in particular of roots) it is by 
default assumed that Nc is at least 3. However, this may be 
manually changed by setting the option NcMin, for example 
CDot[Cl, C2, NcMin -> 100]. 


CDot[{Cl,. . . , Cn}, Options] 
Options with default value: 
NcMin 3 
Verbose — > False 


Similar to CDot [CI, C2] but calculates the scalar product 
of a set of color tensors and themselves and returns the result 
as a List where the ith element is CDot[Ci, Ci]. In addition 
to assuming a minimal value of Nc for simplifications, it 
is also possible to get progress information by setting the 
option Verbose to True. 


CDotMatrix[{Cl,. . . , Cn}, Options] 
Options with default value: 
NcMin 3 
Verbose — > False 


Calculates the scalar product matrix given a List of color 
vectors, i.e., element ij is CDot[Ci, Cj]. If the list contains 
a basis, this function thus returns the scalar product matrix. 


CNorm[Cl, Options] 
Options with default value: 
NcMin -> 3 


Calculates the norm of a color tensors using CDot. For 
simplifying the result, it may be useful to set the option 
NcMin to a large value. 


CNorm[{Cl,. . . ,Cn}, Options] 
Options with default value: 
NcMin -> 3 
Verbose — > False 


Similar to CNorm above, but calculates the norms of a List 
of color vectors {CI,. . . ,Cn}, and returns a List containing 
the norms. Progress information is available by supplying 
the option Verbose— > True. 


CGamma[{Cl,. . . ,Cn}, pl,p2, Options] 

Options with dciault value: 

NcMin -> 3 

Verbose — ^ True 

MakeChecks -» True 

BasisType — > GeneralBasis 


Describes the effect of gluon exchange between partons pi 
and p2, on basis vector Cj as a column vector j in a result- 
ing matrix (technically List of List), i.e., element ij is the 
resulting color structure's i-th component. The BasisType 
option should assume one of the values GeneralBasis, Or- 
thogonalBasis, OrthonormalBasis or TraceBasis. The lat- 
ter is defined to be any basis where the basis vectors consist 
of one product of closed and open quark-lines, not a sum. 
By default it is checked that the color tensor, resulting after 
gluon exchange, is the same when squared directly, as when 
basis decomposed, and then squared. This checks that the 
basis is complete and that the basis decomposition is correct. 
For orthonormal bases, it is also checked that the resulting 
matrix is symmetric. Using MakeChecks — > False these 
checks are turned off. By default, progress information is 
also written out, which can be changed with Verbose — > 
False. 
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